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The cosmic coincidence problem is a serious challenge to dark energy model. We suggest a 



quantitative criteria for judging the severity of the coincidence problem. Applying this criteria to 
three different interacting models, including the interacting quintessence, interacting phantom, and 
interacting Chaplygin gas models, we find that the interacting Chaplygin gas model has a better 
chance to solve the coincidence problem. Quantitatively, we find that the coincidence index C for the 
interacting Chaplygin gas model is smaller than that for the interacting quintessence and phantom 
models by six orders of magnitude. 
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The type la supernovae observations in 1998 discovered that the expansion of the universe is accelerating This 
discovery of late time cosmic acceleration imposes a big challenge and provides opportunities to gravitational and 
particle physics. To explain the accelerated expansion, we can modify either the left hand side (modify the theory of 
gravity) or the right hand side (add an exotic component with negative pressure, dubbed as "dark energy" ) of Einstein 
equation. Although lots of dark energy models have been proposed in the literature, the nature of dark energy is still 
1 a mystery. For a review of dark energy models, see for example [2[ and references therein. 

The simplest candidate of dark energy which is consistent with current observations is the cosmological constant. 
■ Because of the many orders of magnitude discrepancy between the theoretical predication and the observation of 
vacuum energy, the origin of the smallness of the value of the cosmological constant is unknown. Furthermore, the 
. ■ cosmological constant faces the "coincidence problem", namely, why the energy density of dark energy and dark 
matter happens to be of the same order now? In terms of the parameter r, 

o ; 

0^ . A Pdm s 

O. r== ' 1 

Pde 

where p is the energy density, the subscript de denotes dark energy and dm labels dark matter, the coincidence 
problem says why r becomes order of 1 now. If the current value of r which is order of 1 is independent of initial 
conditions, then the coincidence problem is alleviated. Therefore, the resolution of the coincidence problem lies on the 
attractor solution. Along this line of reasoning, the coincidence problem has been extensively studied [H, 0, H, [fil • For 
the quintessence model without interaction, the attractor solution with acceleration is the scaling solution with r = 0, 
so the interaction between dark matter and dark energy is proposed to get nonzero r attractor solution. Assuming 
that the interaction is turned on recently, the standard radiation and matter dominated eras can be recovered and 
the coincidence problem is alleviated because r starts from order 1 unstable fixed point [|| . 

To solve the coincidence problem, r should not vary too much through the whole history of the universe, in addition 
to having attractor solution. Thus, during most of the history of the universe, dark matter and dark energy evolve 
almost in the same way. In other words, during the matter domination, dark energy behaves like dark matter. 
The Chaplygin gas model satisfies this requirement At early times, the Chaplygin gas model behaves like dark 
matter; and it behaves like dark energy at late times. Furthermore, the interacting Chaplygin gas model was shown 
to have the attractor solution with nonzero r in [8j|. So it can be used to solve the coincidence problem. Here the 
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Chaplygin gas model is used as a dark energy model, not as a model unifying dark matter and dark energy. Recently, 
a relationship between the Hubble parameter H and the ratio r was found in [9(. Since the evolution of the Hubble 
parameter H(z) can be determined by observations [l(J EH, the evolution of the ratio r can be derived. In this letter, 
we discuss the evolution of r for several interacting dark energy models. 

The letter is organized as follows. In section II, we study the coincidence problem in the interacting quintessence 
(IQT) model, the interacting phantom (IPT) model and the interacting Chaplygin gas (ICG) model. We conclude 
the letter with some discussions in section III. 



II. THREE INTERACTING MODELS 



Due to the interaction between the two dark components, the equations of motion for pd e and pdm become 

pde + SH'YdePde = -I\ (2) 

and 

Pdm + SHjdmPdm = I\ (3) 

where a dot stands for the derivative with respect to the cosmic time t, the barotropic index 7 is defined as 

7 = 1 + - = 1 + w, (4) 
P 

w is the equation of state parameter, jdm — 1, and T is the interaction term between dark energy and dark matter. 
The origin of the interaction between dark energy and dark matter is not clear and the interaction can be rather 
arbitrary. Here we choose the interaction term [J, |5| 

T = 3Hc(p dm + Pde), (5) 

where the constant c is the coupling constant. A possible mechanism for this interaction from scalar-tensor theory of 
gravity is discussed in detail in [8|, H2, LL3| . 



A. The evolution of r in three interacting models 

Combining Equations @, ([3]) and ([5]), one obtains the variation of r with respect to the cosmic time, 



r = 3Hr 
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The evolution of r for the IQT and IPT models with constant 7 has been discussed in 0, S Q ■ The attractor solution 

is mil 



w 



+ \7w 2 + 4cw 



w — Vw 2 + Acw ' 

with the condition < c < |w|/4. For the ICG model, p c h — —A/p c h, so 

7«* = l-^(l + r) a A', 



(7) 



(8) 



where A' = A(3Hq)~ 2 k 4 and k 2 — 8ttG. The interacting Chaplygin gas model was studied in detail in [8|. The 
attractor solution is r s = c/(l — c) with the condition < c < 1. 

With the help of Friedmann equation, the differential equation of r with respect to H reads Q , 



dr 
dH 
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In terms of the dimensionless variable u, v which are defined as 



U = (3H 2 ) 1 K 2 p dm , 



(10) 
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FIG. 1: r = u/v as a function of ln(l + z) in ACDM, IQT, IPT and ICG models. The initial conditions and the parameter 
are taken as: u| 2= o + Qb = 0.3, v\ z= o = 0.7, and Qb = 0.0389. 
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equations (J5J) and {3]| become 



du 

dx 



-3u + 3c(u + v), 
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dv 
dx 



-3v(l + w de ) - 3c(u + v), 



where x = In a = — ln(l + z). The Friedmann equation becomes 

H 2 

—2 =u + v + n b (i + zf 



(13) 
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where ilb is the fraction of the baryon energy density. The above equations describe the evolution history of the 
universe. We solve the above equations numerically for the ACDM, IQT, IPT, and ICG models, and the results are 
shown in FigureQ] The parameters in the equations are taken as follows. The Hubble constant Hq = 70 kms _1 Mpc _1 
[Til ]. The baryon energy density fib = 0.0389 [H]]. The equation of state parameter w q = —0.9 for the IQT model 
and w p — —1.1 for the IPT model. The coupling constant between dark energy and dark matter c q = 1.0 x 10~ 4 for 
the IQT model, c p = 1.0 x 10~ 4 for the IPT model, and c Cp = 0.04 for the ICG model. Note that for the IQT model, 
observations require c < 2.3 x 10~ 3 [6], so we take take c = 1.0 x 10~ 4 @. 

Note that there is no interaction between the baryon and dark energy because the standard model of particle physics 
is well constrained and the coupling between dark energy and ordinary matter is strongly constrained by the solar 
system experiments [lj|. From figure [T] it is clear that the ICG model solves the coincidence problem because the 
ratio r does not change much during the evolution of the universe and the current ratio r s is the attractor. In the low 
redshift region, the three models almost behave in the same way. But in the high redshift region, the ratio r in the 
ICG model is effectively a constant, but it continually increases to a very high value for the IQT and IPT models. 

Currently there is no data available for the evolution of the ratio r(z) yet. However, by solving equation ([9]), we 
can express H as a function of r, H = H(r). For some special cases, especially for the quintessence and phantom 
models with constant equation of state parameter w, the explicit function forms were obtained in [9J. Recently, the 
Hubble data H(z) become available by using the age of the oldest galaxies PUHH, an( l ^ has been used to constrain 
cosmological parameters [la |. 

In Figure 2, we show the evolution of H with respect to r. Form Figure 2, one sees that in the ICG model, the 
ratio r is almost a constant for H/Hq > 2, while in the IQT and IPT models, r increases with respect to H . 

In figure 3, we show the evolution of H with respect to the redshift z. 
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FIG. 2: r — u/v as a function of ^ in ACDM, IQT, IPT and ICG models. The parameters are taken as the same as that of 
figure 1. 
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FIG. 3: H as a function of z in ACDM, IQT, IPT and ICG modls. The parameters are taken as the same as that of figure 1. 
The 9 data of H(z) (denoted by circles) with la error bars are taken from [hj, while the other 15 data of H(z) (denoted by 
points) with la error bars are taken from [ill ]. 

From Figures 1, 2, and 3, we see that the ICG model is a much better model for overcoming the coincidence 
problem. But in what degree? We need a quantitative criteria for the selection. We shall discuss this point in the 
next subsection. 

B. A quantitative criteria for coincidence 

The interacting models are often invoked to solve the coincidence problem. As we have mentioned before, the 
attractor behavior in the interacting model is a necessary ingredient to overcome the coincidence problem. For the 
attractor solution, the ratio r is independent of the initial conditions, and therefore greatly softens the coincidence 
problem. For such models, the ratio only depends on the parameters of the model, but is independent of the initial 
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conditions. Therefore there is no problem of fine-tuning the initial conditions. However, we still need to choose 
appropriate model parameters to get the right value of r s . In this sense, there exists the problem of fine-tuning the 
parameters. If the ratio r does not change much during the whole history of the universe for most of the parameters, 
then the coincidence problem can be solved. Based on these considerations, we suggest a quantitative criteria for 
judging the severity of the coincidence problem. We divide the coincidence problem into two smaller problems. The 
first is the coincidence between the value of r e at early time and that of ro at present. The second is the coincidence 
between ro and the attractor value (if it exists) rf. To describe the coincidence problem quantitatively, we define the 
indices of early coincidence and late coincidence, respectively, 

C e 4 ^, (15) 



C f 4 r -L. (16) 

The closer to 1 C e or Cf is, the better the coincidence problem is overcome. However we can not define the index 
of the coincidence as C e C/, since C e and Cf maybe varies in the opposite direction, which makes the coincidence 
problem less severe. For example, if C e = 10 10 and Cf = 1CP 10 , the C e C/ = 1. Therefore, we introduce a function, 

F ^ = {l/ X , 0<1<1 W 

By using the function F(x), we can define a proper index of coincidence C in the whole history of the universe as 
follows, 

C ± F{C e )F(C f ). (18) 

which evades the above mentioned problem. 

We must take a standard epoch in the early universe to determine C e . A most natural choice is z — > oo. However, 
we should not expect a phenomenological model can describe the whole evolution history of the universe, especially 
the physics of the very early universe. In this letter, we take z = 100 as "early universe" . 

From equation (|6|), the fixed point of r is 

w + c {rs + 1)2 =0, (19) 

which yields, 



Since r should be real, so 



^ = -d + ^)±V7 + 4?- (20) 



~ > 0, (21) 



or 



2c 



— < -2. (22) 
2c " V ' 



The condition (|2~lj) leads to negative r s , which is not physical. Hence, for the existence of the fixed point, the 
parameters w and c are required to satisfy the condition (|22p . If w is a function of cosmic time, like the case in 
the Chaplygin gas model, then in (|2"2"|) we take the the value of w as lirrij—,—! w. For the ICG model, we take 
c = 0.06, A 1 = 0.4, which is consistent with observations For the IQT model, we take c = 1.0 x 10 -4 Q. 

We apply this criteria to the three interacting models. The result is shown in table I. From table I, it is clear that 
the ICG model evades the coincidence problem. Quantitatively, the coincidence index C for the interacting Chaplygin 
gas model is smaller than that for the interacting quintessence and phantom models by six orders of magnitude. 
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TABLE I: Coincidence indices for the IQT, IPT and ICG models. 



III. CONCLUSION AND DISCUSSION 



In order to study the coincidence problem, we propose a quantitative criteria for the determination of the severity 
of the coincidence problem. 

First, we study the evolution of the ratio r of dark matter to dark energy in the IQT, IPT, and ICG models. Though 
presently we do not have the information about r(z), we can express r with the Hubble parameter H. We found the 
evolutions of r with respect to the redshift z and the Hubble parameter H for the three interacting models. We also 
found the evolution of the Hubble parameter H (z) and compared it with the observational data. We find that the 
ICG model solves the coincidence problem. 

Through a detailed analysis of the coincidence problem, we suggest a quantitative criteria for the determination of 
the severity of the coincidence problem. Applying this criteria to the IQT, IPT, and ICG models, we find that the 
ICG model solves the coincidence problem and the coincidence index C for the interacting Chaplygin gas model is 
smaller than that for the interacting quintessence and phantom models by six orders of magnitude.. 
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